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ABSTRACT 

Planetary systems discovered by the Kepler space telescope exhibit an intriguing feature. While 
the period ratios of adjacent low-mass planets appear largely random, there is a significant excess of 
pairs that lie just wide of resonances and a deficit on the near side. We demonstrate that this feature 
naturally arises when two near-resonant planets interact in the presence of weak dissipation that 
damps eccentricities. The two planets repel each other as orbital energy is lost to heat. This moves 
near-resonant pairs just beyond resonance, by a distance that reflects the integrated dissipation they 
experienced over their lifetimes. We find that the observed distances may be explained by tides if tidal 
dissipation is unexpectedly efficient (tidal quality factor ^ 10). Once the effect of resonant repulsion 
is accounted for, the initial orbits of these low mass planets show little preference for resonances. This 
could constrain their origin. 



1. INTRODUCTION 

NASA's Kepler mission is revolutionizing our knowl- 
edge of planetary systems. It has already discovered 
thousands of transiting planetary candidates, includ- 
ing hundreds of sy stems with two or more planets 
(jBatalha et al.l l2012f ) . Most of these arc Neptune- or 
Earth-sized planets. To date, one of the most intrigu- 
ing Kepler discoveries is that, while the spacing between 
planets appears to be roughly random, there is a dis- 
tinct excess of planetary pairs just wide of certain res- 
onances, and a nearl y empty gap ju st narrow of them 
(jLissauer et all 120111 : iFabrvckv et ani2012[ ). These fea- 
tures are particularly prominent near the 3:2 and 2:1 res- 
onances, and a ffects planets that fall within a few percent 
of resonances ()Fabrvckv et al.|[20Tl . 

Is this resonance asymmetry a feature planetary sys- 
tems are born with, or one they acquire much later on? 
Many studies have reported that planets become trapped 
into first-order reson ances wh en th ey migrate in proto- 
plane t ary disks (e.g..lLee fc P calc 2 00^ iSnellgrove et al.l 
120011: iPapaloizou fc Szuszkiewiczl |2005|) . In fact, the 
presence of resonances among giant planets detected by 
radial velocity has be en regarded as strong evidence 
for d isk migration (e.g.. lMarcv et al.ll2001l : pTinnev et al.l 
l2006f ). However, Kepler's low- mass planets appear to be 
less influenced by resonances, and the pile-ups just out- 
side resonances are partly counterbalanced by the gaps 
inside them. 

In this paper, we identify a process that can modify 
the pair separation and give rise to the observed reso- 
nance asymmetry. But first, let us consider a commonly 
invoked mechanism, tidal circularization. If the inner 
planet is eccentric, tides raised on it would damp its ec- 
centricity, decrease its semi-major a xis, and hence in- 
crease the period ratio of th e pair (jNovak et al.l 120031 : 
iTerguem fc Papaloizoul l2007tFI . Adopting the equilib- 
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rium tide expression from lHu"3 (jl981| ). the damping rate 
for a psudo-synchronized planet is 
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where q = M^/rrii is the mass ratio of the star to 
planet, k2 the tidal love number, i?i the inner planet's 
radius and a its orbital separation. In this tidal model, 
Ti = R\/{GmiTi) where ti is the assumed constant tidal 
lag time which we take to be ti = Pi /(2Qi ), with Qi 
the in ner planet's tidal quality factor (jGoldreich fc Soterl 
Il966( ) and Pi its orbital period. Numerically, 
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The orbital decay rate is a /a = 2ee because orbital 
angular momentum is largely conserved (assuming that 
e ^ 1). So tidal evolution could have potentially circu- 
larized orbits inward of ^ 10 days. As it does so, it moves 
the inner planet inward by 5a /a ~ This increases 

the period ratio for a planet pair by a fractional amount 
of 3ef/2 = 1.5%(ei/0.1)^. However, tidal circularization 
alone can not reproduce the observed asymmetry: as- 
suming all near-resonant pairs were initially uniformly 
distributed in their period ratios, all systems march to 
larger period ratios by a comparable amount. This pro- 
duces neither gap nor peak. 

A more selective mechanism is required. In this paper, 
we show that for a pair of planets that happen to lie near 
a mean-motion resonance, dissipation causes the planets 
to repel each other. The rate of repulsion is greatest at 
exact resonance and falls off steeply away from resonance. 
Planets that are initially slightly closer than resonance 
are pushed wide of the resonance; those that are initially 
wider are pushed even further apart. And planet pairs 
far away from the resonance are not affected. So the com- 
bined action of resonant interaction and damping natu- 
rally give rise to the observed resonance asymmetry. This 
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effect, w hicti we term "res onant repulsion," was invcsti- 
gated bylGreenberd (I1981D for tlie Galilean satellites; by 
iLithwick fc Wul (I2008D to possibly account for the or bits 
of Pluto's minor moons; and by Papaloizouj (|20H|) for 
multiple planet systems. iBatvgin fc Morbidellil ( 20121 ) 
independently arrived at many of the results presented 
in this paper; their paper was posted to arxiv.org at the 
same time as this one. 

2. RESONANT REPULSION 

We consider the evolution of two planets orbiting a star 
and assume that the interaction between the planets is 
predominantly due to the 2:1 resonance. We will also 
include weak external eccentricity-damping forces. The 
energy (or Hamiltonian) of the two planets is, to leading 
order in eccentricity. 



ing terms. 
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(/lei cos (2A2 - Ai - wx) + 7262 cos (2A2 - Ai - 072)) (3) 

where we foll ow ^standard notation (e.g., 

[Murray fc DermottI |2000[ ). with the orbital parame- 
ters for the inner planet denoted by {ai, ei, Ai, roi}, 
and those for the outer planet subscripted by 2. The 
mass of the star and planets are Af*, toi, ''tt-2, and the 
Laplace coefficients are /i = —(2 -I- aD/2)h\i2 
/2 = (3/2 + aD/2)b\^2 - 2a ()Murrav fc DermottI [2OOOI ). 

Near 2:1 resonance (a = 2^^/^), the Laplace coefhcients 
are A = -1.19 and /2 = 0.428. 
We choose units such that 



GM, = 1 



(4) 



and assume that the eccentricities are small. In terms of 
the complex eccentricity 

z, ^ e,e^^^ , (5) 

the equations of mot i on for planet .7 a re (e.g. 
iMurrav fc DermottI [2000HLithwick fc Wull2008D 
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To leading order in mj/M^,, the semi- major axes are 
constant, and the equations for Xj are 
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is the fractional distance to nominal resonance. When 
A < the pair is on the near side of resonance, other- 
wise it is on the far sideQ The 7ej in Equations (fT^ - (fT3|) 
denote the eccentricity damping rates on each of the two 
planets due to some external force (e.g., tides or a dissi- 
pative disk). We assume that ^ |A7t,2|. 

We discard the free solutions to Equations ([T^ - p^ 
because they decay to zero at the rates 7ej, much faster 
than the rate of semi-major axis evolution, as we shall 
see below. The forced eccentricities are, to first order in 
7ej/(An2) < 1: 
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The small phase shift, 0{'yej), relative to the undamped 
forced eccentricities plays a crucial role in resonant re- 
pulsion. 

Inserting the above forced ecc entricities into the sem i- 
major axis equations yields, as in lLithwick fc Wul ()20Q8D . 
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and we have included additional damping terms with 
rates laje'j- This form for the damping rate is appli- 
cable for any process that conserves angular momentum, 
such as tides (see below). By contrast, if the planet is 
migrated in a disk, or pushed by tides raised on the cen- 
tral body (as for Jupiter's moons), the induced rate of 
change of a would be independent of eccentricity. We 
shall not consider those kinds of forces. 

Wc conclude that the distance to resonance changes at 
the rate 



dA 
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The eccentricity equations become, after adding damp- 



* For brevity, we often refer to nominal resonance (A = 0) 
as simply resonance. This should not be confused with a pair 
being locked in resonance, i.e. in a state where the resonant angles 
librate. 
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Fig. 1. — Effect of resonant repulsion on the period distribution 
of planet pairs. The black histogram sh ows Kepler data for planet 
pairs near the 2:1 and 3:2 resonances HBatalha et al.l 12011 1. The 
three red curves show the effect of resonant repulsion on an initially 
flat distribution of pairs at three later times. The parameters are 
chosen such that Amig = .005{t/Gyr)^''^ in Equation I I25I I. similar 
to our fiducial values for tidal damping (Eq. I I26I I). The pileup 



and the evacuated region extends to ' 



where 
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for |A| ^ 1. We verify this rate with an N-body simula- 
tion below. 

As long as r > 0, as we shall argue is the case, then 
A always increases, independent of the sign of A. A 
pair of planets that is initially spaced closer than nom- 
inal resonance (A < 0) will tend to be pushed outside 
of resonance, i.e. to A > 0. And a pair initially out- 
side of resonance will be pushed even further apart. We 
term this effect resonant repulsion. Furthermore, since 
the speed of migration is slowest far from resonance, the 
region near nominal resonance (A = 0) should be unoc- 
cupied, and resonant pairs should evacuate the resonance 
region and pile up outside. This will lead to an asym- 
metry, with more planets outside of nominal resonance 
than inside. 

Two planets that initially have A = Aq repel each 
other to A > Aq, and at time t they migrate to 



where 
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Figure [T] illustrates the effect on the distribution of period 
ratios. 

The sign of F is always positive due to eccentricity 
damping alone, i.e. to the ^ej terms in Equation p3)). 



Furthermore, if tides are the source of damping, then 
-faj — 27ej by angular momentum conservation, leaving 
F > 0; this is also true for any form of damping that 
conserves angular momentum. Other forms of damping 
could in principle result in values of "faj that make F 
negative. However, the fact that Kepler pairs are piled 
up outside of resonances argue that this did not happen. 

3. RESONANT REPULSION BY TIDES 

In this section we focus on the case when the dissipa- 
tion is provided by tidal damping. The rate of eccentric- 
ity damping 7ei is given by Equation ([2]). In addition, 
7ai = 27ei by angular momentum conservation, and we 
may ignore tides on the outer planet (7e2 = 7o2 = 0) be- 
cause tidal damping rates are steep functions of orbital 
period. Therefore Equation becomes 
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Figure [2] shows an N-body simulation with tides of two 
planets initially on the near side of resonance. Resonant 
repulsion pushes them to the far side, in agreement with 
the analytic solution (Equations (|24)) and (|26l) ). There 
is modest disagreement when the pair crosses through 
nominal resonance when the expansion in small e be- 
comes invalid. The free eccentricities damp away after a 
brief initial period (;^ 2 x 10^ yr). On crossing nominal 
resonance, they are regenerated, but then quickly damp 
away again. Damping locks the system into libration 
(of both resonant angles), but this has little dynamical 
significance, as it is merely a consequence of the eccen- 
tricities taking on their purely forced values. 

Figure [3] shows the "resonant repulsion time" (<:rr) 
for all reported Kepler pairs. This is the timescale over 
which resonant repulsion by tides moves a pair towards or 
away from the nearest first order resonance. Mathemat- 
ically, tRR, = |A/A|, where A is the observed fractional 
distance and A is the rate predicted by resonant repul- 
sion (Equation (|22p ) assuming tidal damping is operating 
with Qi ~ 10, k2 = 0.1, and using the observed planet 
and stellar parameters. On this plot, systems that have 
<RR longer than their age have not experienced signifi- 
cant resonant repulsion, while all those with shorter tRR 
should have moved to the right. 

A number of inferences may be drawn. First, most sys- 
tems far from resonances (|A| > 10%) have experienced 
negligible resonant repulsion and were most likely born 
with the period ratio they have today. 

Second, systems within 1 — 10% of resonance exhibit 
fRR that are as long as, or longer than, the typical age of 
systems (a few Gyrs). This is consistent with resonant 
repulsion by tides: systems with shorter ipR would have 
been moved to the right until iRR was comparable to the 
age of the system. Near the 2:1 resonance, it appears 
that pairs as far left as 1.8 and as far right as 2.2 could 
have been affected by the repulsion. 

Last, many systems very near resonances (|A| < 1%) 
exhibit such short tRR that they should have migrated 
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Fig. 2. — An N-body simulation of resonant repulsion, where 
the dissipation is provided by tidal damping on the inner planet. A 
pair of planets initially on the near side of the 2:1 resonance (period 
ratio P2/P1 = 1.99) is pushed to the far side as a result of the inner 
planet moving inward and the outer planet outward. The planets 
both have mass IOMq and orbit a solar mass star, with Pi = 5 
days. To speed up the simulation, we artificially enhance tides by 
assuming a radius of 12Re for the inner planet, while Qi = 10 and 
kj =0.1. The simulation was performed with the SWIFT package 
l|Levison fc DuncanI 1199"^ . modified to include routines for t idal 
damping and relativistic precession. IILevison fc Duncan|[T99ll . 

to much larger A values. At first sight, their presence 
is troubling. However, an inspection of the Kepler cata- 
logue reveals that many of these are in triples or higher 
multiple systems, and these planets are engaged simul- 
taneously in two or more 2-body resonances. The worst- 
off cases are in simultaneous resonances, reminiscent of 
the L aplace resonance of Jupiter's moons (|Yoder fc Peald 
I1981D . Moreover, the fraction of multiples is much higher 
amongst systems with ^rr falling below the solar age line 
than for other random pairs. Our simple picture of reso- 
nant repulsion fails when the planet is subject to two or 
more resonances. In this case, exact resonance may be 
maintained for a much longer time because the planets 
form a heavy ladder with an effectively large inertia. The 
prevalence of simultaneous resonances in these short ^rr 
systems spurs us to hypothesize that all pairs with short 
^RR in Fig. [3] are results of 3-body effects; and that these 
resonances are not primordial, but a combined effect of 
resonant repulsion and 3-body effects. 

Removing the colored circles in Fig. |31 we see a rela- 
tively clear picture that most pairs stay where they were 
born with, while pairs very close to resonances experi- 
ence repulsion and are shifted by a few percent to larger 
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Fig. 3. — The timescale for resonant repulsion to move the period 
ratio of Kepler planet pairs by a distance |A|, where A is the 
observed fractional distance to the closest first order resonance. 
We adopt KIC syst em parameters, with updated values for KOI- 
961 (red dots) from [Muirliead et al. (2012). For tidal dissipation, 
Ql = 10 and k2 = 0.1. The lower panel zooms in to the resonant 
region. If t^iYi ^ system age (the horizontal line is the age of 
the Sun), the period ratios should have evolved little since birth; 
while for tj^R <^ age, we do not expect the systems to linger at 
the observed ratios. The fact that most pairs lie at or above the 
horizontal line is consistent with resonant repulsion by tides. Close 
inspection of systems with very small tRR reveal many are related 
to 3-body effects: the turquoise circles indicate pairs where one or 
both planets are engaged in at least two resonances simultaneously 
(defined as |A| < 3%). Our simple picture of resonant repulsion 
may break down in these cases. 'Uncertain systems' refer to those 
where the nominal total mass > lOOOMgj (assuming Earth density) 
and we discard them from consideration for fear of contamination. 



period ratios. 



4. DISCUSSION 



In this work, we investigate the peculiar fact that there 
is an excess of Kepler planet pairs just wide of resonance, 
and a deficit just inward of resonance. We propose that 
dissipation is responsible for this asymmetry. Two nearly 
resonant planets whose eccentricit i es are weakly damped 
repel each other (lGreenberd ll981l: iLithwick fcWul 120081 : 
iPapaloizoul 120 lit) . This is because dissipation damps 
away the planets' free eccentricities, but the eccentric- 
ities that are forced by the resonance persist despite dis- 
sipation. Planets are typically repelled when dissipation 
acts on these forced eccentricities. As such, resonant 
interaction allows dissipation to continuously extract en- 
ergy from the orbits. Resonant repulsion pushes pairs 
from the near side to the far side of resonance, and nat- 
urally explains the Kepler result. Pairs accumulate at 
a fractional distance Amig wide of each resonance, with 

Amig ^ (A^^i/idamp)^/^, whcrc ^damp is the typical eccen- 
tricity damping time and t the system age (Equations 

(El-®)- 

For the source of dissipation, we focused on tidal damp- 
ing in the inner planet. The typical distance planets can 
repel each other is of order a few percent or less for Ke- 
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pier parameters if the tidal damping is efBcient. The 
deficit of pairs immediately inward of resonance may be 
explained by this repulsion. And the distances outward 
of resonance where planet pairs are found are consistent 
with the theoretically estimated repulsion distance. 

However, a number of inconsistencies between theory 
and data require further investigation. For instance, 
many pairs remain very close to resonance despite a short 
resonant repulsion time. These are often found in sys- 
tems with more than two planets where the planet pairs 
are engaged simultaneously in more than one resonance. 
Wc therefore speculate that in fact all systems with short 
resonant repulsion time are consequences of 3-body ef- 
fects. This may be confirmed using transit-timing varia- 
tion or other tools. 

If resonant repulsion is the reason behind the resonance 
asymmetry, its signature should be observable in future 
studies. The planets should currently have nearly zero 
free eccentricities, and as a result both of the resonant 
angles should be locked at their center-of-resonance val- 
ues, with very small libration amplitude. This can be 
tested with radial velocity measurements or with transit- 
time variations (pthwick et al., ,2012i) . Further more, if 



tidal damping is the dominant dissipation mechanism, 
we expect that the resonance asymmetry should vanish 
for planets at orbital periods greater than 10 — 20 days. 
Long-term Kepler monitoring will decide between tides 
or alternative damping mechanisms, e.g., damping by a 
gaseous or planetesimal disk. 

Our study suggests that the initial period distribution 
of Kepler planets was relatively flat, without major pile- 
ups at or near resonances This is in contrast to jo- 
vian mass planets and could help constrain the origin of 
these low-mass planets. If disk migration is responsible 
for their current location, it must somehow have avoided 
pushing the planets into resonances, perhaps because the 
migration rate was very fast — faster than the resonant 
libra tion rate. Alternatively, planets may be formed in- 
situ ([Hansen fc Murravll2011l ) and have therefore avoided 
convergent migration. 
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^ The chain of resonances observed in systems like KOI-500 , 
KOL730 ijLissauer et al.l [20Tll 'l. KOL2038 ijFabrvckv et al.l [20111 



might also not be primordial, but a combined result of resonant 
repulsion (a 2-body effect) and 3-body interactions. 



